Integer Quantum Hall State in Two-Component Bose Gases 
in a Synthetic Magnetic Field 
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We study two-component (or pseudospin-1/2) Bose gases in a strong synthetic magnetic field. 
Using exact diagonalization, we show that a bosonic analogue of an integer quantum Hall state 
(without intrinsic topological order) appears at the total filling factor ;/ = 1 + 1 when the strengths 
of intracomponent and intercomponent interactions are comparable with each other. This provides a 
prime example of a symmetry-protected topological phase in a controlled setting of quantum gases. 
The real-space entanglement spectrum of this state is found to be comprised of counter-propagating 
chiral modes consistent with the predicted edge theory. 
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There has been an ever growing interest in artificially 
created gauge fields in ultracold atomic gases. While 
atoms are charge neutral, gauge fields can be induced by 
rotating gases iQ or by optically dressing atoms [^, || . A 
sufficiently strong synthetic magnetic field for ultracold 
atoms is expected to offer interesting analogues of quan- 
tum Hall states with a variety of statistics and internal 
states of constituent particles. A relevant parameter for 
describing a strong magnetic field regime is the filling fac- 
tor ly — N/N^, where N is the number of atoms and A^^ is 
the number of magnetic flux c^uanta piercing the system. 
For scalar bosons, a variety of quantum Hall states have 
been predicted at various integer and fractional i^(< 6) 
iQ such as a bosonic Laughlin state ai v ~ 1/2 |j] and 
non-Abelian Read-Rczayi states [^ at i^ = fc/2 with in- 
teger fc> 2 H. 

From a general perspective, Laughlin and Read-Rezayi 
states are examples of topologically ordered states — in 
the sense that they host fractional excitations in the bulk 
and exhibit ground-state degeneracies that depend on 
the system's spatial topology 0. Recently, there has 
been an intensive interest in a new type of phases which 
are not topologically ordered but are distinct from triv- 
ial phases as long as certain symmetries are imposed. 
Such phases are termed symmetry-protected topological 
(SPT) phases, and examples include topological insula- 
tors [Q and Haldane chains |^. Very recently, a system- 
atic classification of SPT phases for interacting bosons 
has been proposed |lO|, ^, opening up possibilities of 
various SPT phases which were not known before. It has 
been predicted that a system of bosons subject to the 
particle-number conservation can exhibit an analogue of 
an integer quantum Hall (IQH) effect whose Hall conduc- 
tivity is quantized in units of an even integer |ll| . It is 
interesting to ask whether such new SPT phases can be 
realized in quantum gases with tunable interactions and 
gauge fields. Quite recently, Senthil and Levin Q have 
proposed that two-component Bose gases in a magnetic 
field may realize one of such phases — a bosonic IQH 
state with its Hall conductivity quantized to 2. This in- 



triguing proposal was based on an effective field theory, 
aiming to discuss a general possibility without specify- 
ing a microscopic Hamiltonian. An important issue of 
the competition among the bosonic IQH state and other 
candidate states in a realistic Hamiltonian has remained 
to be settled. 

In this Letter, we numerically study a simple realis- 
tic system of two-component (or pseudospin-1/2) Bose 
gases in a strong magnetic field. Recent exact diago- 
nalization studies of this system using a torus geometry 
p3| , p^ have identified a series of incompressible states at 
V = fc/3 -|- fc/3 (fc/3 filling for each component) with inte- 
ger fc, in the case of pseudospin- independent interactions. 
Possible candidate quantum Hall states for these filling 
factors include non-Abelian spin-singlet states with an 
SU{3)k symmetry [^; such states have been found to 
appear for fc = 1 and 2 [Q Q. The nature of incom- 
pressible ground states (GS) at higher z/ has remained 
to be identified. Here, by combining calculations on 
torus and sphere geometries, we analyze the competition 
among candidate quantum Hall states. We show that the 
bosonic IQH state of Rcf. 12 does appear at i/ = 1 -|- 1 



when the strengths of intracomponent and intercompo- 
nent interactions are comparable with each other. This 
result establishes the appearance of a SPT phase in a 
simple microscopic model, and opens up a unique possi- 
bility of its realization in a controlled setting of ultracold 
atoms. 

We consider a system of a two-dimensional (2D) Bose 
gas (in the xy plane) having two hyperfine spin states 
(labeled by a =t, i) and subject to a synthetic magnetic 
field B along the z axis. In the case of a rotating system, 
the magnetic field B = 2MVl/q is induced in the rotat- 
ing frame of reference, where M is the particle's mass, 
il is the rotation frequency, and g is a fictitious charge 
for neutral atoms. Unique vortex structures for mod- 
erate rotation frequencies have been investigated both 
theoretically llfl, O and experimentally |lq|. Quantum 
Hall states are expected to emerge as the number of flux 
quanta (or vortices), N^, becomes comparable with the 
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FIG. 1: Candidates for incompressible GSs in the {N^,N) 
plane, calculated on a sphere geometry in the SU{2)- 
symmetric case g^^ — g. Filled circles indicate GSs with the 
total angular momentum L = 0, where incompressible states 
can appear; the area of each symbol is proportional to the 
neutral gap A„. Empty circles indicate the GSs with L > 0. 
Solid and broken lines indicate the relation A'^ = v{N^ + 5) for 
5 = and 2, respectively. Data points are missing for large 
A*'^ and N due to exponentially long computation time. Be- 
low the dotted line with v — 2/3, GSs are highly degenerate 
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number of atoms, iV, for sufficiently strong B. For such 
large B, it is natural to assume that the motion in the 
xy plane is restricted to the lowest Landau level (LLL). 
Within this subspace, we consider the interaction Hamil- 
tonian consisting of intracomponent and intercomponent 
contact interactions: 
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where Na is the number of particles in state a, r"'s are 
the 2D positions of such particles, and g and g^i are 
the effective coupling constants in the 2D plane. We 
mainly focus on the case of 5'f/(2)-symmetric interactions 

.9ti = .9! ^^^ *^^s^ of 9ti 7^ 9 '^il^ '•'Iso ^^ discussed. 

To study bulk properties, it is useful to work on closed 
uniform manifolds having no edge. In Refs. |1^, |lj], a 
periodic rectangular geometry (torus) of sides L^ and Ly 
is used. Here, we also use a sphere geometry and combine 
the two geometries to obtain geometry-independent con- 
clusive results. For a sphere geometry [|l^, a magnetic 
monopole of charge N^{2Trh/q) with integer N^ = 25 
is placed at the center, and it produces a uniform mag- 
netic field B on the sphere of radius R = iV^, where 
i = y/h/{qB) is the magnetic length. Introducing spher- 
ical coordinates {0,ip), single-particle orbits in the LLL 
are given by ipm oc uS+m^s-m ^j^j^ ^ _ sin(0/2)e*'^/^ 
and V = cos(e'/2)e-*'^/2, where m = -5', -S* + 1, . . . , S* is 
the z-component of the angular momentum. On both the 



sphere and the torus, the filling factor in the thermody- 
namic limit is given hy 1/ = N/N^ with N = N^ + N^. On 
a finite sphere, however, the relation between N and iV^ 
involves a characteristic shift S for incompressible states: 
N = v{N^ + (5), where 5 depends on a given candidate 
wave function. Therefore, on a sphere, competing in- 
compressible states leading to the same v in the thermo- 
dynamic limit can be studied separately with different 
{Nci,,N) if they have different shifts. On a torus, there 
is no shift, and all candidates for the same i^ compete in 
the same finite-size calculation. While a torus can give 
less biased results, a sphere is more useful in discussing 
relative stabilities of given candidates. Henceforth, we 
take the units in which h= 1 and £ = 1. 

In Fig. 0, we set (7-1-4^ = g > and look for candi- 
dates of quantum Hall states in the (A^,^, N) plane using 
a sphere geometry. Because of the rotational symme- 
try of a sphere, many-body eigenstates can be classified 
by the total angular momentum L. It is known that in- 
compressible states in general appear as unique GSs with 
L = 0, which are indicated by filled circles in Fig. Q^. The 
area of each filled circle is proportional to the neutral gap 
A„ to be defined later. Solid and broken lines indicate 
the relation N = v{N^ + 5) for (5 = and 2, respectively. 
The bosonic IQH state with v = 2 has 5 = Q |12]; the 
SU{Z)k states with v — 2A:/3 have 6 — 2 and can appear 
only if iV0 + 2 is a multiple of 3 ||l|l. We find that L = 
GSs appear for all [N^^ N) corresponding to the IQH and 
SU{i)k states. Furthermore, we also find L = GSs with 
appreciable energy gaps along the 5 = lines with v ~ A 
and 6, whose nature will be discussed later. Below we 
analyze the stabilities of these candidate states and, in 
particular, the competition between the IQH state and 
5'C/(3)3 state at i' = 2 based on the data of energy gaps. 

We consider two types of energy gaps, charge and 
neutral gaps. The charge gap is defined as Ac{N) = 
Eo iN/2 + 1, N/2) + Eo (iV/2 - l,iV/2)-2£;o (7V/2, N/2) 
for even integer N{> 4), where Eo{N^, N^) is the ground- 
state energy for given particle numbers {Nf,Ni). The 
neutral gap A„ is defined as the gap to the first excited 
state for (iVt'^4.) = iN/2, N/2). 

Figure ||(a) presents the result for a sphere with N^ = 7 
and varying N. Both the IQH state and the SU{3)k 
states can appear in this case. The two types of gaps, 
Ac/2 and A„, show essentially similar behaviors [the di- 
vision by 2 is done in the former because of a doubled 
energy scale involved in the definition of Ac] . The most 
prominent gap is found at iV = 6, where the Halperin 
(221) state ^ [or the SU{3)i state] gives the exact GS 
[ pi[ . The second largest gap (comparable with the first) 
is found at A^ = 14, which corresponds to the 1/ = 2 IQH 
state. This gap is (by a factor sa 1.5) larger than the gap 
at iV = 18, which corresponds to the SU{3)3 state. This 
indicates that of the two candidates at :^ = 2, the IQH 
state is likely to be more stable. It is useful to compare 
Fig. p|(a) with the previous result on a torus [O, pi. 



0.06 



I 0.04 





£ 0.02 



qa. 




f SU{3)2 
v=4/3 



SU(3)3 
v=2 



SU(3)4 
v=8/3 



10 



10 12 14^^16 1? 



20 22 24 



to 

>. 

0) 
LU 



0.1 

,0.08 

0.06 

0.04 

0.02 





.(b)v=2 ^ 


<1 III I> O . 
1 . 1 . 1 . 1 . 


■ (c) v=4 0. 

, , , , 1 , . . 


- ^^^^B-^gf"^ 

. Aj,/2 (sphere) ^^ \^__ 
An (sphere) -a- 
Aj /2 (torus) -B- 

An (torus) -V; 







0.05 y^o.^ 



0.15 .|^,^ 0.05 



FIG. 2: (a) Charge and neutral energy gaps, Ac and A,i, cal- 
culated on a sphere with A^^ — 7. (b,c) The size-dependence 
of the energy gaps at i^ = 2 and 4 on the sphere and torus 
geometries. For a torus, we set L^ = Ly (periodic square). 
For a sphere, we set 5 = 0; namely, calculations are performed 
along the solid lines in Fig. hi The dotted lines indicate ex- 
trapolation of the charge gap to the thermodynamic limit, 
using the data with A'^ > 10. 



where comparable gaps sX v ^ 2/3 and 2 were found. 
Therefore, interpreting the v ^ 2 GS as the IQH state 
gives a better consistency between the sphere and torus 
data. In Fig. ||(a), we also find a large gap at TV = 20 
(comparable with the two largest), whose implication is 
currently unclear. 

In Fig. H(b) and (c), we examine the size-dependence 
of the gaps at i^ = 2 and 4 on the sphere and torus ge- 
ometries. A sphere tends to give larger gaps than a torus 
for finite sizes. For i/ = 2 in Fig. ||(b), we find that the 
difference between the sphere and torus data decreases 
as N increases. Extrapolation to the thermodynamic 
limit (dotted lines) yields Ac/{2g) sa 0.44 and 0.49 for 
the sphere and torus cases; they are in good agreement, 
given the remaining oscillations in finite-size data used 
for the extrapolation. Furthermore, the torus data give 
an additional support for excluding the appearance of 
the 5'C/(3)3 state. If the SU{3)3 state is stabilized in 
the thermodynamic limit, A„ in the current definition 
should drop because of the 10-fold topological GS degen- 
eracy expected on a torus |l^ (a real neutral gap should 
be found above the degenerate GSs in this case). We 
find no appreciable drops of A„ for iV = 6, 12, 18 (sys- 
tem sizes where the SU{3)3 state is allowed) in Fig. 0(b). 
For z^ = 4 in Fig. 0(c), the system size is not sufficient 
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FIG. 3: (Color online) Real-space entanglement spectra of 
the i^ = 2 IQH state with (A^^, N) = (8, 16) on a sphere. We 
plot the entanglement energies {— logpi} defined from the 
eigenvalues {pi} of pA, where A is the northern hemisphere; 
the lower the entanglement energy is, the larger weight it has 
in pA- The entanglement energies are classified by the number 
of t and 4- atoms, A'^^ and N^^ , and the z-component of the 
total angular momentum, L^ , on A. 



to check the consistency between the torus and sphere 
data; yet, simple linear extrapolation of the charge gap 
(dotted lines) indicates a tendency toward a finite value 
in the thermodynamic limit on both geometries. 

For 1/ = 2, we have also analyzed the case of g-^i ^ g 
by calculating the spectra on a torus |22]. We find that 
the IQH state with an appreciable energy gap survives 
against moderate changes in the ratio g-fi/g- 

We have so far discussed the appearance of the IQH 
state at 1/ = 2 using only spectral properties. A useful 
way to unveil the topological features of this state is to 
investigate the entanglement spectrum p3[. This is an 
eigenvalue spectrum of a reduced density matrix pA on 
a subsystem A, and has been demonstrated to serve as a 
useful probe of the edge spectrum. Here we use the GS 
on a sphere, split it into two hemispheres, and calculate 
the real-space entanglement spectrum (RSES) Q for the 
northern hemisphere A. The result is presented in Fig. ^. 
Below we analyze this result on the basis of an effective 
field theory. 

The i^ = 2 bosonic IQH state is described by the 
Chcrn-Simons theory with the 2 x 2 K matrix with 
Ki2 = K2i = l and K^ = K22 = Q |l|. Through a 
standard procedure Q] , the effective Lagrangian for one- 
dimensional edge modes is obtained as 

-C = --—{Kal3dt(t)adx<j>l3 + K,/? 9j; 0a c'x 0/3 ) , (2) 

47r 

where ■^dx4'a describes the density of bosons of spin 
a =t, \, and Vap is the velocity matrix. Diagonalizing C 
via charge and spin modes (f)^/g = (0-|-±0j_)/v2, one finds 
that the charge mode is right-moving while the spin mode 



is left-moving |Q. Introducing the velocities v^jg > for 
these modes, the Hamiltonian and the total momentum 
for the edge of length L^ are given by 



H = 



P 



^[vc{d.cl>cf+vs{d,cl>,f], (3a) 



dx 
An 



[id,(bcf - {d^sf 



(3b) 



By performing mode expansions for the bosonic fields 
4)c/si the terms in the above equations can be diagonal- 
ized as 



j-{dAc,s) =-r 



Att 



Lx 



(AiV^ ± A7VJ 



T7l — 1 



''rn,c/s 



(4) 



Here, ANa {a =t: i) is a change in the number of spin-a 
bosons relative to the GS; {um.c} and {rim^s} are sets of 
non-negative integers describing oscillator modes, and re- 
spectively give rise to the well-known f/(l) boson count- 
ing of degeneracies: 1,1,2,3,.... 

In Fig. H, the lowest entanglement energy is found at 
Lf = -16 with {n\,N\) = (4,4), which can be iden- 
tified with the "GS" of virtual edge modes probed by 
the RSES. "Excited states" in the RSES are found in 
both positive and negative directions of Lf relative to 
the "GS", which indicates the counter-propagating na- 
ture of the edge modes. This is in contrast to the case 
of a single chiral edge mode [p^ pil, where low- lying 
"excited states" are found only in one direction of Lf. 
The excitations along the left envelope of the spectrum 
can be interpreted as the left-moving spin mode with 
P < 0, as spin excitations with {AN\, AN\) = (+1, -1) 
and (+2, —2) relative to the "GS" are found there; the 
shifts AL^ = — 1 and —4 of these excitations are con- 
sistent with Eqs. ^ and (^. Similarly, the excitations 
along the right envelope can be interpreted as the right- 
moving charge mode with P > 0, as charge excitations 
with (AnI^AN^) = (-1,-1) and (-2,-2) are found 
there. The levels appearing between the two modes can 
be identified with combinations of these modes. 

Finally, we discuss the nature of the gapped states at 
1/ = 4, 6, . . . found in Fig. |. The RSES of the z/ = 4 state 
exhibits a counter-propagating nature of edge modes [^ , 
as in the v = 2 case in Fig. b[ This suggests that the v = 4 
state has certain similarities with the v = 2 IQH state. 
However, the bosonic IQH states with v > A predicted in 
Ref. |lli are not good candidate for the present system, 
as their corresponding simplest K matrices are not sym- 
metric with respect to the interchange of the two com- 
ponents. A more detailed characterization of the states 
at z/ = 4, 6, . . . is left for future studies. 

In summary, we have numerically demonstrated that a 
bosonic IQH state without intrinsic topological order ap- 
pears in a simple realistic model of two-component Bose 



gases in a synthetic magnetic field. Given the recent suc- 
cess of creating a time- reversal- invariant gauge field [p5l , 
it is interesting to further investigate how the symme- 
try protection gives rise to novel topological phases and 
edge states in ultracold atomic systems in a variety of 
synthetic gauge fields. 
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Supplementary Material 
Stability of the IQH state at ;/ = 2 

Here we consider the case of g^i ^ g and discuss 
the stability of the IQH state found at :/ = 2. Fig- 
ure H presents the spectra for {N^, N) = (6, 12) on a 
torus geometry as a function of the coupling ratio .g-i-4,/5 
or its inverse. We find that an appreciable gap above 
a non-degenerate ground state remains in the regions 
0-7 < .9ti/ff — 1 (I'^ft panel) and 1 > g/pti ^ 0.3 (right 
panel). This indicates that the IQH state is stable against 
moderate changes in the ratio g^i/g around the SU(2)- 
symmetric case. 

Entanglement spectra 

Figure || presents the real-space entanglement spectra 
(RSES) of candidate incompressible states on a sphere ge- 
ometry. Here we divide the sphere into two hemispheres, 



calculate the reduced density matrix pA on the north- 
ern hemisphere A, and plot the entanglement energies 
{—log Pi} defined from the eigenvalues {p,} of pA- The 
case of {N^,N) = (6,24) in Fig. ^(a) corresponds to the 
incompressible state with i^ = 4 and (5 = 0, whose na- 
ture has been elusive. Just like the 1/ = 2 case in Fig. 3 
of this Letter, the low-lying entanglement levels appear 
in both positive and negative directions of L^, indicat- 
ing the counter-propagating nature of edge modes. This 
suggests that the v ~ A state has certain similarities with 
the v = 2 IQH state. For comparison, Fig. @(b) presents 
the RSES of the Halperin (221) state [or SU{3)i state] 
at ly = 2/3. The low-lying levels appear only in the neg- 
ative direction of Lf, indicating the chiral (east-moving) 
nature of edge modes of this state. 
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FIG. 4: Energy spectra versus coupling ratio g-fi/g or its 
inverse for {N^,N) = (6,12). Diamonds (O) indicate two 
lowest-energy states at JC = in the equal-population case 
A'^ = Ni = N/2; the ground state energy of this sector is 
subtracted from the whole spectrum. Circles (O) indicate 
other eigenstates in the equal-population case. Crosses (x) 
show eigenstates for the minimally imbalanced case A'^ — 
N/2 + 1 and Aj, = A/2 - 1. 



FIG. 5: Real-space entanglement spectra for {N,/,,N) = 
(6,24) and (10,8) on the sphere geometry. We plot the en- 
tanglement energies {— logpi} defined from the eigenvalues 
{pi} of PA, where A is the northern hemisphere. The en- 
tanglement energies are classified by the numbers of ^ and ^ 
atoms, A-|. and Aj^ , and the jz-component of the total angular 
momentum, L- , on the northern hemisphere. 



